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We present a detailed analysis of the incident-photon-energy and polarization dependences 
of the resonant inelastic x-ray scattering (RIXS) spectra at the Cu K edge in La2Cu04. Our 
analysis is based on the formula developed by Nomura and Igarashi, which describes the spec- 
tra by a product of an incident-photon-dependent factor and a density-density correlation 
function for 3d states. We calculate the former factor using the 4p density of states from an 
ab initio band structure calculation and the latter using a multiorbital tight-binding model 
within the Hartree-Fock approximation and the random phase approximation. We obtain 
spectra with rich structures in the energy-loss range 2-5 eV, which vary with varying momen- 
tum and incident-photon energy, in semi-quantitative agreement with recent experiments. 
We clarify the origin of such changes as a combined effect of the incident-photon-dependent 
factor and the density-density correlation function. 
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1. Introduction 

Taking advantage of strong synchrotron light sources, the resonant inelastic x-ray scatter- 
ing (RIXS) has recently become a powerful tool to probe charge excitations in solids. K-edge 
resonances are widely used in transition- metal compounds, because it could detect momentum 
dependence of charge excitations. 1-16 It is described as a second-order optical process that a 
ls-core electron excites to an empty 4p state with absorbing a photon, then charge excitations 
are created in the 3d states to screen the core-hole potential, and finally the photoexcited Ap 
electron recombines with the ls-core hole with emitting a photon. In the end, the charge 
excitations are left with carrying the energy and the momentum transferred from photon. 

Intensive studies have been carried out on Cu oxide perovskite compounds for better un- 
derstanding of the unconventional high-Tc superconductivity. 3-5,7-16 In the undoped material 
La2Cu04, the spectra are found to be composed of several peaks which changes with changing 
momentum transfer. 5 ' 13 Such spectra have been theoretically studied by several groups. 17-22 
Nomura and Igarashi (NI) 19,20 have proposed a general formalism of the RIXS spectra by 
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extending the resonant Raman theory developed by Nozieres and Abrahams 23 on the basis of 
the many-body formalism of Keldysh. The NI formula makes it possible to calculate the RIXS 
spectra on complicated models including many orbitals, and provides clear physical interpre- 
tations to the RIXS spectra. It is composed of two factors, one describing an incident-photon 
dependence and the other the density-density correlation function for the 3d states. Similar 
formulas which are proportional to the density-density correlation function have been derived 
by using different methods. 24 

This formula has been applied to a two-dimensional cuprate La2Cu04 20 and quasi-one- 
dimensional cuprates SrCu03 19 and CuGe03. 9 In these studies, the electronic structures have 
been calculated on the d-p model within the Hartree-Fock approximation (HFA) in the an- 
tiferromagnetic (AFM) phase. It is known that the HFA works well to describe electronic 
structures in the AFM insulators. Two-particle correlations have been taken into account 
within the random phase approximation (RPA). The calculated spectra as a function of en- 
ergy loss have reproduced well the RIXS spectra varying with varying momentum. The RIXS 
spectra have been interpreted as a band-to-band transition. The RPA correlation had to be 
included for obtaining better agreement with the experiment. 

Subsequently, the present authors have analyzed multiple-scattering contributions due 
to the core-hole potential, because the core-hole potential is not definitely weak. 22 Having 
evaluated the multiple-scattering contributions by means of the time-representation method 
by Nozieres and De Dominicis, 25 we have found that main contributions could be absorbed 
into the shift of the core-level energy with minor effects on the RIXS spectral shape. This 
result partly justifies the use of the Born approximation. Quite recently, we have demonstrated 
the usefulness of the NI formula by analyzing the RIXS spectra in NiO. 26 Using the HFA and 
the RPA on the multiorbital tight-binding model, we have obtained the RIXS spectra, which 
vary with varying momentum, in quantitative agreement with the experiment. 27 

Recently, Lu et al. 13 have presented comprehensive RIXS data with systematically chang- 
ing the incident-photon energy in La2Cu04. Kim et al. 14 have recounted the momentum 
transfer dependence and incident photon polarization dependence of the RIXS line shape in a 
number of cuprates. Ellis et al. 15 have reported the high resolution RIXS spectra in La2Cu04. 
Their data show fine structures as a function of energy loss, in addition to the peaks found 
in previous data, which vary with varying the incident-photon energy. Although our previous 
analysis using the d-p model has been successful to clarify the origin of the RIXS spectra, 20 
the analyses have been limited to the situation of the incident-photon energy corresponding 
to the excitation of the Is electron to the peak of the 4p density of states (DOS). Also, it 
seems hard to explain the newly observed fine structures. Encouraged by the success in the 
application to NiO, 26 we extend the d-p model to a multiorbital tight-binding model, which 
includes all O 2p and Cu 3c? orbitals as well as the full intra-atomic Coulomb interaction 
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Fig. 1. Sketch of the assumed crystal and magnetic structure of La2Cu04 with omitting La sites. 
The tiny tilting, rotation and the Jahn- Teller distortion of CuOg octahedron are neglected. 

between 3d orbitals, and analyze the incident-photon-energy and polarization dependences of 
the RIXS spectra in comparison with the experiment. 13 ' 15 

In the present study, we obtain the AFM insulating solution having an energy gap ~ 1.7 eV 
within the HFA. Note that the band structure calculation with the local density approximation 
(LDA) fails to reproduce the insulating state. These energy bands by the HFA are used to 
calculate the density-density correlation function. We treat the two-particle correlations within 
the RPA. The incident-photon-dependent factor is calculated by using the 4p DOS obtained 
from the ab-initio band structure calculation. We obtain richer structures in RIXS spectra in 
a range of energy loss 2 ~ 5eV than those obtained by the d-p model analysis, 20 which vary 
with varying momentum transfer and incident-photon energy. The obtained spectra are in 
semi-quantitative agreement with the experiments. 13-15 We clarify the origin of such spectral 
change as a combined effect of the incident-photon-dependent factor and the density-density 
correlation function. 

The present paper is organized as follows. In Sec. 2.1, we introduce the multiorbital tight- 
binding model. In Sec. 2.2 we discuss the electronic structure within the HFA in the AFM 
phase of La2Cu04. In Sec. 3, we briefly summarize the NI formula for the RIXS spectra. In 
Sec. 4, we present the calculated RIXS spectra in comparison with the experiment. The last 
section is devoted to the concluding remarks. 
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2. Electronic Structure of La2Cu04 

2.1 Multiorbital tight binding model 

For explaining the fine structures in the observed RIXS spectra, we need to extend a simple 
d-p model to a multiorbital tight-binding model, in which most of parameters are determined 
from an ab-initio band structure calculation. 

As shown in Fig. 1, we approximate the crystal structure of La2Cu04 by a simple tetrag- 
onal one without the Jahn- Teller distortion and the tilt of the Cu06 octahedron. Extending 
the d-p model, we introduce a tight-binding model involving all Cu 3d orbitals and O 2p 
orbitals at the in-plane and apical sites. We exclude orbitals belonging to the La atoms, since 
those orbitals play minor roles in the electronic states near the insulating gap. 28 Thereby the 
Hamiltonian may be written as 

H = H + H h (1) 

^0 = E ima n ima + E P jma n P jrn(J + (timjm'dimaPjm'a + H.C.^ 

irna jrna amm' 

+ EE (ffm,y m >P)rncPj'rn>a + K. C )+Y, E (fii^mAma^'m' a + H.C.) , (2) 

(J J') (7mm ' (M'} amm' 

Hl = ^EE E E E 9{mia 1 m 2 a2]m 3 a 3 m4a 4 )dl mia A im2a2 d im4a J im3a3 . (3) 

i m\a\ m^ai 7713173 m^a± 

The part Hq represents the kinetic energy, where di ma and Pj ma denote the annihilation 
operators of an electron with spin a in the 3d orbital m at Cu site i and that of an electron 
with spin a in the 2p orbital m at the O site j, respectively. Number operators nf mcr and n p ma 
are defined as nf ma = d\ ma d ima , n p jma = p) ma p jm a- The transfer integrals, t^ Jm „ t pp m j , m ,, and 
tfmi'm' are evaluated from the Slater-Koster (SK) two-center integrals, (pda), (pdir), (ppa), 
(ppir), (dda), (ddn), (dd5). 29 We use the SK parameters determined by the ab-initio band 
structure calculation. 28 The part Hi represents the intra-atomic Coulomb interaction on Cu 
sites. The interaction matrix element g {y\Vi\ ^3^4) (v stands for spin-orbit (ma)) is written in 
terms of the Slater integrals F°, F 2 , and F 4 . Among them, F 2 and F A , which are known to be 
slightly screened by solid-state effects, are taken from the cluster model exact diagonalization 
analysis of the X-ray photoemission spectroscopy. 30 On the other hand, F° is known to be 
considerably screened, so that we regard the value as an adjustable parameter. The Coulomb 
interaction on O sites is absorbed into a renormalization of the O 2p level parameters E p ma . 
The d-level position relative to the p-levels is given by the charge-transfer energy A defined as 
A = E Cg — E p + 9U in the d 9 configuration. 31 Here U is the multiplet-averaged d-d Coulomb 
interaction given by U = F° — (2/63) F 2 — (2/63) F 4 . We treat the charge transfer energy A as 
an adjustable parameter in our calculation. We fix the charge transfer energy at A = 2.76 eV 
in order that the edge position of the RIXS intensity is consistent with the experiment. Note 
that the energy gap is determined by the charge transfer energy, since the electron system 
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Table I. Parameter values for the tight-binding model of La2Cu04 in units of eV. SK parameters are 
taken from ref. 28. E £g — E t2g is the crystal field splitting between the Cu 3d e g and ti g states. 
^O(i) an d ^0(2) are t ne O 2p levels for the in-plane and apical oxygen, respectively. Slater Integral 
F 2 and F 4 arc taken from ref. 30. 



Cu 


F° 


7.3, 11.3 


F 2 


11.41 


F 4 


7.308 




E e 9 ~ E t2g 


1.18 


A 


2.71 








ddo 


0.065 


ddir 


-0.067 


dd5 


-0.079 


Cu-O(l) 


pda 


1.253 


pdir 


0.853 






Cu-0(2) 


pda 


0.568 


pdn 


0.377 






0(1) 


ppa 


0.586 


ppir 


-0.384 


^O(l) 





0(l)-0(2) 


ppa 


-0.207 


ppir 


-0.196 






0(2) 


ppa 


0.171 


ppir 


-0.025 




0.50 



in La2Cu04 belongs to the charge transfer type insulator. We have checked that the RIXS 
spectra do not crucially depend on the values of F°. The parameters used in the calculation 
are listed in Table I. 
2.2 Hartree-Fock approximation 

Assuming the AFM ordering, we choose the unit cell that contains two Cu sites, one of 
which is a up spin site and another is a down spin site, and eight O sites, four of which are the 
in-plane 0(1) sites and the others apical 0(2) sites. Then we solve the Schrodinger equation 
H^> = by disregarding the fluctuation terms in Hj, i.e., 

H ? F = \Y1 E r(0) (W> ^4) <4*<W4*<U> ( 4 ) 

i U1U2U3U4 

where is the antisymmetric vertex function defined by 

T (0) (y x v 2 \ ^3^4) = 9 {v x vi\ u 3 u 4 ) - g (y x v 2 \ v^) , (5) 

with the bracket (• • • ) indicating the ground state average. 

For parameter values of F°, we take up two typical values F° = 7.3 and 11.3 eV, which 
correspond to the values of U used in refs. 32 and 20, respectively. For both cases, we obtain 
stable AFM solutions which have the energy gap E gSuP ~ 1.7 eV and the spin moment at Cu 
site scu ~ 0.3 h. The fact that the energy gap depends little on the F° values indicates that 
the system belongs to the charge transfer insulating phase. The density of states and the 
dispersion curves along some symmetric lines are shown in Figs. 2 and 3, respectively. The 
weight of the minority spin x 2 — y 2 states is indicated by the length of vertical bars in fig. 3. 
Labels assigned to the DOS peaks and to the dispersion curves will be related to the peaks 
in the RIXS intensities. 
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Fig. 2. (Color online) Density of states projected on the Cu 3d x 2 — y 2 , 3z 2 — r 2 , and 0(1) 2p a 
states. The origin of energy is at the top of valence band, (a) and (b) are for F° — 7.3 and 11.3 
eV, respectively. 



The states in the conduction band (assigned A) consist mainly of the minority spin x 2 — y 2 
state. Thereby the RIXS intensity is strongly correlated to the minority spin x 2 — y 2 state in 
the valence band, as discussed in the next section. The states around the top of valence band 
(assigned B) are mainly composed of the in-plane O pa states and the majority spin x 2 — y 2 
states, but they also contain small weights of the minority spin x 2 — y 2 states, as seen from the 
figure. Their weights strongly depend on the momentum k. Particularly, in the valence bands 
located just below the gap, the weight is negligible along the (tt/2, tt/2)-(tt, 0) line, while it 
becomes discernible along the (0, 0)-(7r/2, 7t/2) and (tt, 0)-(0, 0) lines as denoted by B. The 
states labeled C and D noticeably contain the minority spin x 2 — y 2 states. We note that the 
states relating to the DOS peaks and the dispersion curves denoted by A, B, C, and D do not 
strongly depend on the F° values. The energy position and the weight of the x 2 — y 2 states are 
mainly controlled by the charge transfer energy A and the SK parameters. We also note that 
while the dispersion curves for the states denoted A and B are quite consistent with those 
calculated by the d-p model, it is difficult to find the correspondence for the other dispersion 
curves between the multiorbital model and the d-p model. 

3. Formula for RIXS 

We briefly summarize the NI formula for the RIXS, following ref. 26. In the RIXS process, 
the incident photon is absorbed by exciting a Cu Is core electron to the unoccupied Cu 4p 
state, and a photon is emitted by recombining the Ap electron and the core hole. This process 
may be described by 

H x = w -i= £ e^p'l^c^e^ + H.c, (6) 

q« V 1 irjij 
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Fig. 3. (Color online) Dispersion curves along the several symmetric lines. The weight of minority 
spin x 2 — y 2 states is shown by the length of the vertical bar. The origin of energy is at the top of 
valence band, (a) and (b) are for F° — 7.3 and 11.3 eV, respectively. 

where w represents the dipole transition matrix element between the Is and the Ap states. 
We assume that w is constant, since it changes little in the energy range of 20 eV above 
the absorption edge. The e^ 1 represents the 77-th component (77 = x, y, z) of two kinds of 
polarization vectors (a = 1,2) of photon. Annihilation operators p\ and Sj CT are for states 
4p v and state Is at Cu site i, respectively. The annihilation operator c qa is for photon with 
momentum q and polarization a. In the intermediate state of the RIXS process, the core-hole 
potential is acting on the 3d states, creating an electron-hole pair. The interaction is described 
by 

Hl s -3d = V Yl 4m* d im*4v> S i*'' ( 7 ) 
imua 1 

where i runs over Cu sites. In the end of the process, an electron-hole pair is left with mo- 
mentum and energy q = (q, u) = (q« — q/,u;j — loj), where qi = (qi,uJi) and qj = (qj,ujf) are 
energy-momentum of incident and scattered photons. 

The RIXS intensity is derived on the basis of the Keldysh-Green function scheme. The 
diagrammatical representation of the process is shown in Fig. 4 (see ref.26 for details). Within 
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Fig. 4. Diagram for the RIXS intensity within the Born approximation for the Is core-hole potential. 
The wavy lines represent photon Green's functions. The solid lines with the labels 4p and Is 
represent the bare Green's functions for the 4p electron and the Is core electron, respectively. The 
dotted lines represent the core-hole potential V. The shaded part represents the density-density 
correlation function of the Keldysh-type. 



the Born approximation to the core-hole potential, we obtain 



W(q u a>i;qf,a>f) 



N\w\ q 



(8) 



Here N is the number of unit cell, (ef v ) is the r\ component of the polarization vector e« 
(e/) with 7] = x,y, z. 

The factor | Y2 V e^L 7 ^ (ujf, to) ef v \ 2 describes the incident-photon dependence, which is 
given by 

V [~> Pi; (e) d6 



(9) 



, eo + e ls + iT ls -e)(oJi-uj + e ls + iT ls - e) ' 
The Vis represents the life-time broadening width of the core-hole state, and eo indicates the 
energy at the bottom of the Ap band. The pj£ is the DOS matrix in the p symmetric states, 
which may be given by 

P% (0 = E E tiri ( k ) *rM ( k ) 6 ( £ " ^ ( k )) » ( 10 ) 

where 4> v j (k) is the amplitude of the Cu p v component in the band state specified by the 
band index j and momentum k with eigenenergy ej (k). The off-diagonal components (rj ^ r/) 
are negligible because of symmetry. This expression comes from the upper triangle in Fig.4. 

The factor Y" d + ~ (q) in Eq. (8) is the density-density correlation function of the Keldysh 
type, which is defined by 



r«T-(q,") = E E *a£,aw(9)> 



(ii) 



\m<T X'm'a' 
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where 

/oo 
<(p q A' m vOVW mCT (0))e^ dr, (12) 
-oo 



PqAma - J "77 ^ d k+qAm CT rf kAm CT • (13) 



with 

*_,__= J 

■ " k 

Here k runs over the magnetic first BZ. The index Amu specifies a tight-binding orbital at 
site A with orbital m and spin a. We assign A = 1, 2 to two Cu sites in the unit cell. 
Using the solution within the HFA, we obtain 

Y + ' {0) (a) ~ n + ~ (0) (a) (14) 

I \ma,\'m'a' W ~ ll \rno\mcr,\' m' o' \'rn' a' W ' V 1 ^ 

with 

n + ~ (0) (a) 

\im\cri\2m2U2,\' 1 m' l u' 1 \' 2 m' 2 cr' 2 

= ^EE*( w - E H k+( i) + E i ( k )) t 1 - n / ( k + <*)] ^ ( k ) 

X <P\ im inj' ( k + q) ^Aimiffij' ( k + ^ ^A' 2 m^4,i (k) ^A 2 m 2 a 2 ,j ( k ) ■ (15) 

where Ej (k) and rij (k) are the eigenenergy and the occupation number of the eigenstate 
specified by jk, respectively. The <p\ ma ,j (k) represents the amplitude of the 3d orbital and 
spin ma at the Cu site A in the energy eigenstate specified by jk. This expression describes 
the excitation from occupied states to unoccupied states with Cu 3d amplitudes. 

We treat the correlation effect on the electron-hole pair by the RPA. It is included into a 
vertex function: 

W = E E A 66,6^)n^ (9) A^ U < (9) • (16) 
66 66 

with the vertex 

A^(q)=\l-tF-(q)]~ 1 (17) 

To simplify the notation, we abbreviate the indices Xma as £. In Eq. (17), / represents a unit 
matrix, and T is the bare four-point vertex given by 

f = T (0) (m 1 a 1 m 4 a4;m2(J2m 3 a 3 )5x 1 x 2 Sx s xJx 1 x 3 , (18) 

■ J 66:66 

with £i = Aimi<7i, £2 = A2W.2O2, £3 = A3?n,30"3, £4 = Xiju^o^ (Aj = 1,2). Note that T is 
non-zero only for £1, £2, £4 belonging to the same Cu site. The two-particle propagator 
F (q) is given by 

(9)1 f e , , = TF E ^4m 4( T4,j (k) ^A 2 m 2 a 2 ,j ( k ) fXimmJ' 0* + O) ^msos J' ( k + ^ 



6-6,66 A , 

k 



(k) [l - rij/ (k + q)] n,-, (k + q) [1 - rij (k)] 



a; - (k + q) + Ej (k) + i«5 u- Ey (k + q) + Ej (k) - id 



(19) 
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Fig. 5. The LDA Ap DOS convoluted by the Lorentzian function with FWHM 2Pi s = 1.6 eV. In 
connection with the absorption coefficient, the origin of energy is shifted so that the prominent 
peak for the p z DOS locates at 8991 eV. A thin line represents the absorption coefficient with 
taking account of the dipole matrix element evaluated by the band calculation. 



For the details of the derivation, see refs.22 and 26. The RPA correction has been found to 
have an important role in the momentum dependence of the spectra. 

4. Calculated Results 

In order to calculate the incident-photon-dependent factor \L^{ojf, uj)\ 2 , we need the Ap 
DOS, pH(e). We evaluate the Ap DOS in a nonmagnetic (NM) phase on the basis of the 
muffin-tin KKR band structure calculation method within the LDA. It is known that the 
LDA fails to predict an insulating phase in the NM phase. However, the 3d states forming a 
metallic phase hybridize little with the Ap bands, giving rise to a minor effect on the Ap DOS 
distributing around 5-30 eV above the Fermi level. Figure 5 shows the Ap DOS convoluted 
with a Lorentzian function of FWHM 2Fi s = 1.6 eV. Reflecting the layered structure, pfp(e) 
is quite different from p%p(e). Note that the Ap DOS is nearly proportional to the Cu -fT-edge 
absorption spectra under the condition that the dipole matrix element is constant and the 
interaction is neglected between the core hole and the Ap electron. In connection with the 
absorption coefficient, we set the energy difference between the Cu Is level and the prominent 
peak in the Ap z DOS to be 8991 eV. 

Figure 6 shows the contour plot of |Lg((jj;t<j)| 2 calculated by substituting the Ap DOS into 
Eq. (9). The enhancement starts earlier with increasing incident-photon energies in the out-of- 
plane polarization (77 = z), reflecting the difference between pfp(e) and (e). Except for these 
points, the general tendency of the enhancement factor looks similar in both polarizations; a 
large peak appears first in the small energy loss region and moves toward the higher energy loss 
region with increasing incident photon energy. The overall intensities decrease with increasing 
incident photon energy. 

Another factor Y^~{q) mainly determines the RIXS spectra as a function of energy loss. 
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Energy Loss CO (eV) 

Fig. 6. Contour plot of \L% (u>i, u>)\ 2 as a function of energy loss w and incident photon energy u>i in 
the out-of-plane polarization (77 — z) (a) and in the in-plane polarization (77 = x) (b) with T\ s = 1 
eV. Units of the intensity arc arbitrary. 

We calculate this factor from Eqs. (14) and (16) with replacing the 5 function by a Lorentzian 
function (FWHM = 0.2 eV) in order to simulate the instrumental resolution. Figure 7 shows 
the calculated results for q along symmetry lines. We show the results with F° = 11.3 eV 
(corresponding to U = 10.7eV in the d-p model analysis 20 ). We have checked that another 
choice F° = 7.3 eV leads to essentially the same result, except for the appearance of a weak 
intensity around ui = 6eV. We obtain continuous spectra ranging from ui = 2 eV to 6 eV. 
Intensities around u = 2, 3.2, and 4.5 eV are caused by charge excitations of B— >A, C^A, 
and D— >A transitions, respectively, within the x 2 — y 2 symmetry in the minority spin states 
(see Figs. 2 and 3). 

We have a prominent peak around 4.5 eV, which stays at the same position with changing 
momentum. We also find that the spectral shape in the low energy region changes with 
changing momentum; a broad hump exists around uj = 2 — 4 eV at q = (0,0), which is 
enhanced by the RPA correction. This hump grows up to become a peak around 3.2 eV with 
changing q along symmetry lines (0, 0) — (tt, 0) and (0, 0) — (tt/2,it/2). Meanwhile, the 3.2 eV 
peak is suppressed with changing q along a symmetry line (tt/2,tt/2) — (vr,7r), vanishing at 
q = (vr,7r). The RPA correction help suppress further the intensity in a region of to = 2 — 4 
eV at q = (n, it), in contrast to the enhancement of the intensity at q = (0, 0). 

The RIXS spectra are given by the product of |L^(u;.;, uj) | 2 and (q). Figure 8 shows 
the calculated spectra as a function of energy loss uj in the out-of-plane polarization, in 
comparison with the experiments. 13 ' 15 At q = (0,0), a hump-like intensity, which is found 
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Fig. 7. The correlation function Y/ (q, uj) as a function of energy loss uj. Top panels are for q along 
(0, 0) — (tt, 0), while bottom panels are for q along (0, 0) — (w, tt). 



around u = 2 - 4 eV in ^ + ~(q = 0, uj), is enhanced by \L z B (uJi,uj)\ 2 with u>i = 8990 - 8994 
eV. Since the enhancement is larger with lower uj, intensities around uj = 2 eV become a 
peak-like shape in the RIXS spectra. Intensities around uj = 3 — 4 eV in Yj (q = 0,ui) are 
moderately enhanced by \L z B (ui, uj)\ 2 , forming a very broad peak in the RIXS spectra. This 
broad spectral feature, which is not produced by the d-p model analysis, 20 is consistent with 
the experiments. 13-15 The large 4.5 eV peak in F d + ~(q = 0, uj) loses its weight in comparison 
with the intensities around w = 2 — 4 eV because of small enhancement by \L z B (u>i;uj)\ 2 . With 
increasing Wj(> 8994 eV), the RIXS intensities around ui = 2 — 4 eV decrease, since the peak 
in \L z B (u>i;uj)\ 2 as a function of energy loss uj moves to a region of higher u. The intensity of 
the 4.5 eV peak first changes little and then decreases with further increasing uii, since the 
\L z B (uJi, lu)\ 2 becomes smaller. With changing q along (0,0) — (tt, 0), and (0,0) — (tt/2,tt/2), a 
peak develops around uj = 3.2 eV in accordance with the development of the peak in Yj (q, uj) 
for u>i = 8992 — 8994 eV. The intensity of the peak around ui = 3.2 eV becomes comparable to 
the 4.5 eV peak at q = (tt, 0) and (tt/2, it/2), while intensities in the region of uj < 3 eV become 
smaller in comparison with the corresponding spectra at q = (0, 0). This behavior corresponds 
well to the experimental line shape at q = (tt,0), which looks like two peaks around uj = 3.2 
and 4.5 eV. With changing q along (tt/2, tt/2) — (tt,tt), intensities around ui = 2 — 4 eV are 
suppressed by the RPA correction in accordance with the change of Y d +_ (q, uj), and only one 
peak is overwhelmingly left at 4.5 eV for a wide range of uii = 8990 — 8996 eV. This explains 
the experimental spectra at q = (tt,tt), which looks like a single peak around uj = 4.5 eV. 
Although the spectral shape is more consistent with the experiments 13-15 than that calculated 
with the d-p model on the whole, several discrepancies between the present calculation and the 
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experiments still remain; the intensity around to = 2 eV is remarkably enhanced at q = (0, 0) 
with tJi = 8990.5 eV, the 4.5 eV peak is hardly discernible at q = (0, 0), the intensity at to = 4 
eV is enough large forming a peak-like structure at q = (0, 0) with uji = 8992 eV, and peak-like 
structures are found in the higher energy loss region oj > 5 eV in the experiments. We note 
that peak structures are observed around ui = 4.5 eV at q ^ (0, 0), which seem to correspond 
well to the calculated 4.5 eV peaks. The calculated spectra at q = (0, 0) specifically fail to 
reproduce the observed spectral structure in u = 4 — 4.5 eV. To remove these discrepancies, 
we may need to take account of the electron correlations beyond the RPA and the effects 
beyond the Born approximation to the core-hole potential. 

Figure 9 shows the calculated RIXS spectra as a function of energy loss u in the in-plane 
polarization. The scale of the calculated intensity is the same as that in the out-of-plane 
polarization. The enhancement by \L B (u>i; lo)\ 2 starts to be effective at u>i about 3 eV higher 
with magnitude much smaller compared to that by \L B (uJi; lo)\ 2 . Reflecting these differences, 
the RIXS spectra are enhanced for uJi > 8993 eV with smaller intensities than those in the 
out-of-plane polarization. Except for these points, the change in spectra with changing q 
and the incident-photon energy are similar to those in the out-of-plane polarization. This is 
consistent with the recent experiment by Kim et al. 14 

In comparison with our previous studies based on a d-p model with fixing uji at 8991 eV 
in the out-of-plane polarization, 22 the spectral change obtained in the present paper with 
changing q is consistent with our previous study. Difference is that the present study using a 
multiorbital tight-binding model leads to a distribution of intensities around 2 — 4 eV. This 
corresponds well to the experimental spectral shape. The present study also improves peak 
position from our previous value u> = 5.5 eV 20 ' 22 to 4.5 eV, although the spectral shape still 
needs to be improved. 

5. Concluding Remarks 

We have analyzed the incident-photon-energy and polarization dependences on the RIXS 
spectra in La2Cu04 on the basis of the formula developed by Nomura and Igarashi. This 
formula expresses the RIXS spectra by a product of the density-density correlation function 
and the incident-photon dependent factor. To explain the fine structures found in the RIXS 
experiment, we have extended the d-p model used in our previous analyses to a multiorbital 
tight-binding model, which includes all the Cu 3d and O 2p orbitals as well as the full Coulomb 
interaction between 3d orbitals. We have calculated the density-density correlation function 
using the HFA and RPA, and the incident-photon-dependent factor using the 4p DOS from 
the ab initio band structure calculation. Fine structures are hardly reproduced by the d-p 
model. In the present detailed analysis, they are found in energy loss 2 ~ 5 eV, which vary 
with varying momentum and incident-photon energy, in semi-quantitative agreement with the 
experiment. Note that the peak shifts as a function of momentum should not be interpreted 
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Fig. 8. (a) RIXS spectra as a function of energy loss in the out-of-plane (z) polarization. Top panels 
are for q along a symmetric line (0, 0) — (ir, 0), while bottom panels are for q along (0, 0) — (n, tt). 
The solid and dashed curves are the spectra calculated within the RP and HF approximations, 
respectively, (b) RIXS spectra reproduced from refs. 13, 15. The curves denoted Ellis are taken 
from the ref. 15. The numbers in panels indicate the incident photon energy u,. 



as a dispersion relation of a kind of exciton, because the peaks with broad widths constitute 
an energy continuum generated by a band-to-band transition. The experimental structures 
are not clear enough for more detailed comparison with the calculated spectra, we hope the 
correspondence will be clarified by improving the instrumental resolution. 

The electron correlation works to modify the single-particle energy bands given by the 
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Fig. 9. RIXS spectra as a function of energy loss for the in-plane (xy) polarization. Each panel is for 
a typical q value. The solid and dashed curves are the spectra calculated within the RP and HF 
approximations, respectively. The unit of intensity is the same as in Fig. 8. 



HFA. 33 One prominent effect is a creation of "satellite" peak around 6 ~ 11 eV below the top 
of the valence band. According to the three-body scattering theory by the present authors, 34 
the creation of satellite works to push the 3d states toward upper energy region in the valence 
band. Such modifications are not strong on the minority spin x 2 — y 2 states in the shallow 
energy region from the top of the valence band. Another effect is the reduction of the energy 
gap from the HFA value. We have taken account of this effect by adjusting the value of A 
to give the experimental energy gap within the HFA. Since the RIXS spectra for u> < 6 eV 
come from the 3d states in the shallow energy region of the valence band, these points may 
explain why the HFA works rather well. As regards the Born approximation to the core-hole 
potential, our calculation explains semi-quantitatively the incident-photon energy dependence, 
indicating that the Born approximation works rather well. Although the present analysis is 
in semi-quantitative agreement with the experiment, several drawbacks remain as pointed 
out in the preceding section. Taking account of electron correlations and going beyond the 
Born approximation to the core-hole potential are desired for quantitative agreement with the 
experiment. 13-15 

Although experimental data have been accumulated for doped cuprates, 10-12 theoretical 
analyses are limited on a one-band Hubbard model within the exact diagonalization method, 35 
and on a three-band Hubbard model analysis within the HFA on the basis of the present 
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formalism. 36 An analysis with a detailed model like the present paper may be necessary 
to clarify the momentum and incident-photon dependences of the spectra. Since electron 
correlations are expected to be more important in doped cuprates, such studies seem rather 
hard and are left in future. 
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